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Abstract 



(^ ■ We recall the main properties of the classical action of Euler-Lagrange 5'ci(x, t;xo), which links 

^ , the initial position xq and its position x at time t, and of the Hamilton-Jacobi action, which connects 

a family of particles of initial action S'o(x) to their various positions x at time t. 

Cn ■ Mathematically, the Euler-Lagrange action can be considered as the elementary solution of the 

Hamilton-Jacobi equation in a new branch of nonlinear mathematics, the Minplus analysis. Phys- 
ic', 
Oj, ically, we show that, contrary to the Euler-Lagrange action, the Hamilton-Jacobi action satisfies 

j-j ' the principle of least action. It is a clear answer on the interpretation of this principle. Finally, 

C^ ■ 

^ ■ we use the relation-ship between the Hamilton-Jacobi and Euler-Lagrange actions to study the 

^\ 

convergence of quantum mechanics, when the Planck constant tends to 0, for a particular class of 

^ I quantum systems, the statistical semiclassical case. 
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I. INTRODUCTION 

In 1744, Pierre-Louis Moreau de Maupertuis (1698-1759) introduced the action and the 
principle of least action into classical mechanics:- "iVainre, in the production of its effects, 
does so always by the simplest means [...] the path it takes is the one by which the quantity of 
action is the least," and in 1746, he states:^ " This is the principle of least action, a principle 
so wise and so worthy of the supreme Being, and intrinsic to all natural phenomena [...] 
When a change occurs in Nature, the quantity of action necessary for change is the smallest 
possible. The quantity of action is the product of the mass of the body times its velocity and 
the distance it moves." Maupertuis understood that, under certain conditions, Newton's 
equations are equivalent to the fact that a quantity, which he calls the action, is minimal. 
Euler,'" Lagrange,^ Hamilton,- Jacobi^ and others, will make this principle of least action 
the most powerful tool to discover the laws of nature.-*^ It allows, with the same approach, 
to determine both the equations of particle motion (if we minimize on the trajectories) and 
the laws of nature (if we minimize on the parameters defining the fields). 

However, for the trajectories of particles, this principle is a problem for many scientists as 
recalled by Henri Poincare, who was nonetheless one of its major users:- " The very statement 
of the principle of least action has something shocking to the mind. To go from one point to 
another one, a material molecule, taken away from the action of any force, but constrained 
to move on a surface, will follow the geodesic line, i.e. the shortest path. It seems that 
this molecule knows the point one wants to lead it to, that it anticipates the time needed to 
reach it along such or such path, and then chooses the most convenient path. In a sense, the 
statement presents this molecule as a free animated being. It is clear that it would be better 
to replace it by a less shocking statement where, as philosophers would say, the final causes 
would not appear to replace the efficient ones." 

We will see that the difficulties of interpretation of the action come from the existence 
of two actions corresponding to two different boundary conditions: the classical action (or 
Euler-Lagrange action) S'ci(x, t; xq), which links the initial position xq and its position x at 
time t, and the Hamilton- Jacobi action ^(x, t), which links a family of particles of initial 
action 5*0 (x) to their various positions x at time t. In the Euler-Lagrange case, the initial 
velocity is unknown and in the Hamilton- Jacobi case, the initial position is unknown. 

In section 2, we recall the main properties of Euler-Lagrange and Hamilton- Jacobi ac- 



tions. In section 3, we see that the Euler-Lagrange action can be considered as the elemen- 
tary solution of the Hamilton- Jacobi equation in a new branch of nonlinear mathematics, 
the Minplus analysis. In section 4, we show that the Hamilton- Jacobi action satisfies the 
principle of least action, contrary to the Euler-Lagrange action, giving a response to the 
scientists' embarrassment concerning the interpretation of this principle. Finally, in section 
5, we use the relation between the Hamilton- Jacobi and Euler-Lagrange actions to study the 
convergence of quantum mechanics, when the Planck constant tends to 0, for a particular 
class of quantum systems, the statistical semiclassical case. 

II. THE EULER-LAGRANGE AND HAMILTON-JACOBI ACTIONS 

Let us consider a system evolving from the position xq at initial time to the position x 
at time t; let x(s) and u(s) be its position and its velocity at each time s G [0,t]. We have: 

^^ = u(s) forsG[0,t] (1) 

x(0) = xo, x(t) = X. (2) 

If L(x, X, t) is the Lagrangian of the system, the Euler-Lagrange action functional is defined 
by 

^£l(u(.)) = / L(x(s),u(s),s)rfs (3) 

where the evolution of x(s) depends on u(s) through equations ([T]) ([2]). 

When the two positions xq and x are given, the Euler-Lagrange action S'ci(x, t;xo) is 
the function which realizes the minimum (or more generally an extremuni^) of the Euler- 
Lagrange action on the velocity field u(.) for all trajectories from (xo,0) to (x,t): 

S'ci(x,t;xo) = minj£;i(u(.)) = min <^ / L(x(s), u(s), s)rfs L (4) 

u{.) u(s),0<s<t Uq J 

the minimum of (jlj) is taken on the continuous controls u(s), s G [0,t], with the state x(s) 
given by the equations ([I])(l2]). 

The solution (u(s), x(s)) of (^, if the Lagrangian L(x, x, t) is twice differentiable, satisfies 
the Euler-Lagrange equations on the interval [0,t]:— 

-^(x(.), x(.), .) - |-(x(s), x(.), .) = (0 < . < t) (5) 

x(0) = Xo, x(t) = X. (6) 



For a nonrelativistic particle in a linear potential field with the Lagrangian L(x, x, t) = 
imx^ + K.x, the equation ([5]) yields ■^{mx{s)) — K = 0. The trajectory which minimizes 
the action is x(s) = xq + f (x — xq) — ^ts + ^s^, and the Euler-Lagrange action is equal 
to ^,Kx, t; xo) = m(^ + ^^^^t - £^t\ 

More generally, in the case of a nonrelativistic particle with the Lagrangian L(x, v,t) = 
|?7iv^ — V(x, t), the Euler-Lagrange action yields the velocities of the two extremities of the 
trajectory: 

Vo = x(o) = --^(x,t;xo) and v(x,t; Xq) = x(t) = — ^(x,t;xo). (7) 

TTl CXq TTl C'X 

For the above exemple, we find: 

x-xo Kt x-xo Xt 
Vo = x(o) = — — and v(x, t; Xq) = x(t) = — h — . 

Let us consider now that an initial action 5*0 (x) is given, then the Hamilton- Jacobi func- 
tional from Xq is defined by:— 

JhAO) = ^o(xo) + / L{x{s), u(s), s) ds = 5'o(xo) + Jbl(u(.)) (8) 

Jo 

where the evolution of x(s) depends on u(s) through equations ([T]) ([2]). 

The Hamilton- Jacobi action S'(x, t) is the function which realizes the minimum of the 

Hamilton- Jacobi action functional on the trajectories arriving in x at time t: 

S'(x, t)= mm Jhj{u{.)) = min < Sq (xq) + / L(x(s), u(s), s)(is > (9) 

xo;u(.) xo;u(s),0<s<t [^ Jo J 

the minimum of (^ is taken on all initial positions Xq, on the controls u(s), s G [0,t], with 
the state x(s) given by the equations ([I])®. 

Because the term 5*0 (xq) has no effect in (Q for the minimization on u(s), we deduce the 
following relation between the Hamilton- Jacobi action and Euler-Lagrange action: 

5'(x, t) = mm{So (xq) + ^^(x, t; xq)). (10) 

xo 

For a particle in a linear potential V{x.) = —K.x with the initial action 5*0 (x) = -mvo ■ x, 
the Hamilton- Jacobi action is equal to S (x, t) = ni\Q ■ x — |?7iVQt -|- K.xt — ^K.vot^ — ^^■ 
The Hamilton- Jacobi action ^(x, t) defined by ([9]) can be decomposed into 

r pt—dt pt 

S{x,t) = min < 5*0 (xo) + / L{x{s),u{s), s)ds + / L(x(s), u(s), s)(is 

xo;u(s),0<s<i (^ Jq Jj_^j 



and then satisfies the optimahty equation: 

S{x,t) = niin < ^(x — / u{s)ds,t — dt) + / L(x(s), u(s), s)(is 

u(s),t-dt<S<t [ J^_^^ _/(_^( 

If we assume S to be differentiable for x and t, L differentiable for x, u and t, and u(s) 
continuous, this equation becomes: 

f dS OS 

= min <^ --- (x,t) u (t) dt - —- (x,t) dt + L(x, u(t),t)dt + o (dt) 
u{t) [ ox at 

and in dividing by dt and letting dt tend towards 0"*", 

— (x,t) = mm <^ L(x,u,t) - ""^ (x,t) > . 

We recall that at all convex functions /(u) : u G M" — > M we can associate its Fenchel- 
Legendre transform /(r) : r G M" — > M defined by /(r) = max(r ■ u — /(u)). The 
Hamiltonian H{x, p, t) is then the Fenchel-Legendre transform of the Lagrangian L(x, u, t) 
for the variable u. 

Then, the Hamilton- J acohi action satisfies the Hamilton- Jacobi equations: 

^ + i/(x,-,)=0 (11) 

S(x,0) = S„(x). (12) 

For the Lagrangian L(x, x, t) = |?7ix^ — l^(x, t), iJ(x, VS*, t) = max(v.VS' — ^ttiv^ + 



\^(x, t)), we have for the optimum rrav = VS. We deduce if(x, VS) = ^ (VS*) + \^(x, t) 

and the well known result: 

The velocity of a nonrelativistic classical particle in a potential field is given for each point 

(x,t) by: 

V x,t = ^^ 13 

m 

where S (x,t) is the Hamilton- Jacobi action, a solution to the Hamilton- Jacobi equations: 
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_ + _(V5)^-H^(x,«)=0 (14) 

S(x,0)=5„(x). (15) 

Equation (TT3l) shows that the solution 5* (x,t) to the Hamilton- Jacobi equations yields the 
velocity field for each point (x, t) from the velocity field — '^^ at initial time. In particular, 
if at initial time, we know the initial position Xinit of a particle, its velocity at this time is 



equal to — °^^'"" . From the solution S* (x,t) to the Hamilton- Jacob! equations, we deduce 
with (fT3|) the trajectories of the particle. The Hamilton- Jacobi action S (x,t) is then a field 
which "pilots" the particle. 

Let us recall how the Hamilton- Jacobi action allows to find Newton's second law 
of motion.— First, we take the gradient of the Hamilton- Jacobi equation (fT^ j^^^ -|- 
- E TT^^ + f^ = 0. Second, we remark that ^ ( ^) = J^ + V j^v^ and with 

m ^-^ axiOx-j ax-j axi ' at \ axi I ataxi ^-^ axiox^ J 

the equation (IT3|) where mv = VS", we conclude ^ (mv) = — VV". 



III. THE EULER-LAGRANGE AND HAMILTON-JACOBI ACTIONS IN MIN- 
PLUS ANALYSIS 

There is a new branch of mathematics, the Minplus analysis, which studies nonlinear 
problems through a linear approach, cf. Maslo'vs^'i^ and Gondran^ii^. The idea is to sub- 
stitute the usual scalar product J^ f{x)g{x)dx by the Minplus scalar product: 

{f,g)= inf {/(a;)+^(x)} (16) 

x£X 

In the scalar product we replace the field of the real number (M, -|-, x) with the algebraic 
structure Minplus (M U {+ oo},min, -|-), i.e. the set of real numbers (with the element 
infinity {-|-oo}) endowed with the operation Min (minimum of two reals), which remplaces 
the usual addition, and with the operation + (sum of two reals) , which remplaces the usual 
multiplication. The element {+00} corresponds to the neutral element for the operation 
Min, Min({-Fcx)}, a) = a\faER. 

This approach bears a close similarity to the theory of distributions for the nonlinear case; 
here, the operator is "linear" and continuous with respect to the Minplus structure, though 
nonlinear with respect to the classical structure (M, +, x). In this Minplus structure, the 
Hamilton- Jacobi equation is linear, because if 5*1 (x, t) and 5*2 (x, t) are solutions to ( TT4|) . then 
min{A -|- S'i(x, t),yU -|- S'2(x, t)} is also solution to the Hamilton- Jacobi equation ( IT4l) . 

The analogue to the Dirac distribution 6{x) in Minplus analysis is the nonlinear distri- 
bution (5min(x) = {0 ?/ X = 0, -|-oo if uot} . With this nonlinear Dirac distribution, we can 
define elementary solutions as in classical distribution theory. In particular, we have: 

The classical Euler- Lagrange action Sci{x,t; ocq) is the elementary solution to the 



Hamilton- J acobi equations l[Tl\)l[T^) in the Minplus analysis with the initial condition 

^(x, 0) = (5min(x-xo) = {0 if X = xq, + cx) ifuot}. 
The Hamilton- Jacobi action ^(x, t) is then given by the Minplus integral: 

5'(x, t) = inf{S'o(xo) + Sc^x, t; xq)} 
in analogy with the solution of the heat transfer equation given by the classical integral: 



^ {x-xpf 

S{x,t) = So{xo)—^e 4:t dxQ. 
J ly-Rt 



In this Minplus analysis, the Legendre-Fenchel transform is the analogue to the Fourier 
transform. This transform is known to have many applications in Physics: this is the one 
which sets the correspondence between the Lagrangian and the Hamiltonian of a physi- 
cal system; which sets the correspondence between microscopic and macroscopic models; 
which is also at the basis of multifractal analysis relevant to modeling turbulence in fluid 
mechanics.— 

IV. THE PRINCIPLE OF LEAST ACTION 

Equation ([9]) shows that, among the trajectories which can reach (x, t) from a position 
at the initial time whose initial velocity field is known. Nature chooses the velocity which 
minimizes (or realizes the extremum) of the Hamilton- Jacobi functional. Then, the principle 
of least action defines the velocity field at time t: (v(x, t) = — ^^^). For the Lagrangian 
L(x, X, t) = imx^ — K.x, we find v(x, t) = Vq + Kt/m. The Hamilton-Jacobi action ^(x, t) 
does not solve only a given problem with a single initial condition ( Xq, — ^^^ ), but a set 
of problems with an infinity of initial conditions, all the couples ( y, — ^^ j . 

Does the Euler-Lagrange action correspond to the principle of least action? The answer 
seems positive from equation (j4]). But in fact, in the absence of an initial velocity field as 
in the Hamilton-Jacobi action, the Euler-Lagrange action answers a problem posed by the 
observer, and not by Nature: "If we see that a particle in Xq at the initial time arrives in 
X at time t, what was its initial velocity vq?" To solve this problem, the observer must 
solve the Euler-Lagrange equations (15|6|) which are difficult because they concern the entire 



trajectory. We have seen that this velocity is given by equation ([7]): vq = — ^-g^(x, t; xq). 
This is an a posteriori point of view. 

Contrary to the Euler-Lagrange action, the Hamilton- Jacobi action answers the following 
problem: "If we know the action (or the velocity field) at the initial time, can we determine 
the action (or the velocity field) at each later time?" This problem is solved from time to 
time by the evolution equation 01 11) or 0141) which is local. This is an a priori point of view. 
It is the problem solved by Nature with the principle of least action. The Euler-Lagrange 
action, which is an elementary solution to the Hamilton- Jacobi equation, seems to satisfy 
this principle. But, only the Hamilton- Jacobi action satisfies the principle of least action. 
It is a clear answer to the Poincare question on the interpretation of this principle. 

V. LIMIT OF THE SCHRODINGER EQUATION IN THE STATISTICAL SEMI- 
CLASSICAL CASE 

Equation (fTOj) between the Hamilton- Jacobi and Euler-Lagrange actions will allow us to 
study the convergence of quantum mechanics, when the Planck constant tends to 0, for a 
particular class of quantum systems. 

Let us consider the wave function solution to the Schrodinger equation \I'(x, t): 

<I'(x,0) = <I'i,(x). (18) 

With the variable change ^^(x, t) = ^yp^'^'x^t) exp(i — |^), the quantum density p''(x, t) and 
the quantum action ^^'(x, t) depend on the parameter h. The Schrodinger equation can be 
decomposed into Madelung equations^ (1926): 

^^'^-'^) + ^iVS\.,t)r + y(x,t) - ^^^S ^ (19) 



dt 2m' ' ' '' ' ' ' 2m ^ p\x,t) 

dp^{x,t) ,, , ., ,V^'^(x,t), ^ ^, , 

L ' +div{p\x,t) ^^^) = V(x,t) 

at m 



(20) 



with initial conditions 

p'^(x,0) =p[;(x) and S\x,0) = S^{x). (21) 

Here, we study the convergence of the density p''(x, t) and the action S'^(x, t), when h tends 
to 0, for a particular preparation of the particles. A quantum system is prepared in the 

8 



statistical semi-classical case^^^ if its initial probability density Po(x) and its initial action 
Sq{x) are regular functions po(x) and 5*0 (x) not depending on h, and its interaction with the 
potential field V(x, t) can be described classically. It is the case of a set of particles that are 
non-interacting and prepared in the same way: a free particles beam in a linear potential, 
an electronic or Cgo beam in the Young's slits diffraction, an atomic beam in the Stern and 
Gerlach experiment. Then, we have the following result:— i^S 

For particles in the statistical semi- classical case, the probability density p^{x,t) and the 
action S^\x,t), solutions to the Madelung equations [W^) [2D\] [2l\) . converge, when h —^ 0, 
to the classical density p{x,t) and the classical action S{x,t), solutions to the statistical 
Hamilton- J acobi equations: 

^ + ci..(p(x,t)^^)=0 V(x,t) (23) 

p(x,0) = po(x) and 5(x, 0) = 5o(x). (24) 

We will demonstrate in the case where the wave function \E'(x, t) at time t is written as 
a function of the initial wave function \l'o(x) by the Feynman paths integral formula3^ (p. 

58): 

*(x,t) = / F(t,n)exp(^5ez(x,t;xo)^o(xo)rfxo 

where F(t, K) is an independent function of x and of xq and where ^^(x, t; xq) is the clas- 
sical action. In the statistical semi-classical case, the wave function is written \E'(x, t) = 



F(t, h) J a/po(xo) exp(|(S'o(xo) + 5'ci(x, t; Xo))(ixo. The theorem of the stationary phased 
shows that, if h tends towards 0, we have \I'(x, t) ~ exp(|?7imx(,(S'o(xo) + S'c/(x, t; Xq)), that 
is to say that the quantum action S'''(x, t) converges to the function 

S'(x,t) = mmxo(S'o(xo) + S'ci(x,t; Xq)) 

which is the solution to the Hamilton- Jacobi equation f l22|) with the initial condition (12^ . 
Moreover, as the quantum density p^{'K,t) verifies the continuity equation (pOl) . we deduce, 
since S^{'x,t) tends towards ^(x, t), that p^{x,t) converges to the classical density p(x, t), 
which satisfies the continuity equation (1231) . We obtain both announced convergences. 

If we consider the system with initial conditions Po(x) = po(x) = (27ro"Q) ae ^"^o and 
5'q(x) = 5'o(x) = mvo ■ x in a linear potential field ^(x) = — K • x, where ao, vq, Co and 



K are constants independent of h, the density p^{'x,t) and the action S^{'x,t) are equal t 

,2 \2 



21 



P (x,t) = (27r(T2(t)) 2exp-^ .,. "'^ and 



c^nr ^ 3^ 1,, . 2^ 1 2 .. 1.. 2 K^^S (x - Co " Vgi - K|^) k^ 



with 0";^ (t) = (To ( 1 + {ht/2mal) j . When /i — > 0, cr^ (t) converges to uq and the den- 

sity p^(x, t) and the action ^^(x, t) converge to p{^.,t) = (27rcrQ)~2e ^''o and 

^(x, t) = — ImVpt + ttivq ■ X + K ■ xt — |K ■ Vot^ — ^|^ which are solutions to statistical 
Hamilton- Jacobi equations (1221) (l23l) (l24l) . 

The statistical Hamilton- Jacobi equations correspond to a set of independent classical 
particles, in a potential field V^(x, t), and for which we only know at the initial time the 
probability density po (x) and the velocity v(x) = — °^^' ' . These particles are not indistin- 
guishable because, if their initial positions are known, their trajectories will also be known. 
Nevertheless, these particles will have the same properties as the indistinguishable ones. 
Thus, if the initial density po (x) is given, and one randomly chooses A^ particles, the N! 
permutations are strictly equivalent and do not correspond to the same configuration as for 
indistinguishable particles. 

For particles prepared in the statistical semi-classical case, the uncertainity about the po- 
sition of a quantum particle corresponds to an uncertainity about the position of a classical 
particle, whose initial density alone has been defined. In classical mechanics, this uncer- 
tainity is removed by giving the initial position of the particle. It would not be logical not 
to do the same in quantum mechanics. It is then possible to assume that for the statistical 
semi- classical case, a quantum particle is not well described by its wave function. One needs 
therefore to add its initial position and it becomes natural to introduce the so-called de 
Broglie-Bohm trajectories^^i^i with the velocity v''(x, t) = — VS'''(x, t). 
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